New precision mass measurements of neutron-rich calcium and potassium isotopes 

and three-nucleon forces 
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We present precision Penning-trap mass measurements of neutron-rich calcium and potassium 
isotopes in the vicinity of neutron number N = 32. Using the TITAN system the mass of J K was 
measured for the first time, and the precision of the 51,52 Ca mass values were improved significantly. 
The new mass values show a dramatic increase of the binding energy compared to those reported 
in the atomic mass evaluation. In particular, 52 Ca is more bound by 1.74 MeV, and the behavior 
with neutron number deviates substantially from the tabulated values. An increased binding was 
predicted recently based on calculations that include three-nucleon (3N) forces. We present a com- 
parison to improved calculations, which agree remarkably with the evolution of masses with neutron 
number, making neutron-rich calcium isotopes an exciting region to probe 3N forces at neutron-rich 
extremes. 



The neutron-rich calcium isotopes present a key region 
for understanding shell structures and their evolution to 
the neutron dripline. This includes the standard doubly- 
magic 48 Ca at N = 28 and new shell closures at N = 32 
and possibly at N = 34. The calcium chain, with only 
valence neutrons, probes features of nuclear forces simi- 
lar to the oxygen isotopes, which have been under inten- 
sive experimental and theoretical studies [TJ [2]. While 
the oxygen isotopes have been explored even beyond the 
neutron dripline, 52 Ca is the most neutron-rich nucleus 
where mass measurements and 7-ray spectroscopy have 
been done. Due to the extremely low production yields 
for N > 32, only the neutron-rich, mid-shell titanium 
and chromium isotopes have been reached to N = 34. 
This makes the semi-magic calcium chain in the vicin- 
ity of N — 32 an important stepping stone towards the 
neutron dripline. 

Phenomenological forces in the pf shell, such as the 
KB3G [3] and GXPF1 g] interactions, have been fit 
to N ^ 32 in the calcium isotopes, but they disagree 
markedly in their prediction for 54 Ca and for a possible 
shell closure at N — 34. In addition, it is well known 
that calculations based only on two-nuclcon (NN) forces 
do not reproduce 48 Ca as a doubly-magic nucleus when 
neutrons are added to a 40 Ca core [5]. Motivated by these 
deficiencies, the impact of 3N forces was recently investi- 
gated in the oxygen and calcium isotopes within the shell 



model [H [7] • It was shown that chiral 3N forces pro- 
vide repulsive contributions to valence neutron-neutron 
interactions. As neutrons are added, these are key for 
shell structure and spectroscopy, and for the determina- 
tion of the neutron dripline. For calcium, the N = 28 
magic number was reproduced successfully, with a high 
2 + excitation energy, and in the vicinity of N = 32, the 
calculations based on NN and 3N forces generally pre- 
dicted an increase in the binding energy compared to 
experimental masses [6]. 

Nuclear masses provide important information about 
the interplay of strong interactions within a nucleus and 
the resulting nuclear structure effects they elicit. The 
systematic study of masses has, for example, led to the 
discovery of the new magic number N = 16 and has 
been key for understanding the island of inversion [8"HTU] . 
Currently, Penning traps |llj provide the most precise 
mass spectrometers for stable and unstable nuclei and 
are in use at almost all rare isotope beam facilities [12] . 

In this Letter, we present the first precision mass mea- 
surements of 51 ' 52 Ca and 51 K performed at TRIUMF's 
Ion Trap for Atomic and Nuclear Science (TITAN) [131- 
IT5] . This presents the first direct mass measurements to 
N = 32 in these nuclides, extending the measurements of 
a previous campaign |16j . For the calcium isotopes, we 
compare the evolution of the new TITAN mass values 
with neutron number to microscopic calculations based 
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FIG. 1. (Color online) TOF-ICR resonances for 51 K, 51 Ca, 
and 52 Ca with excitation times of ~ 80 ms. The applied 
radio-frequency urf is given as the difference from the center 
frequency v c&nt of the scan range. The solid lines are theoret- 
ical fits to the line shape [18| . 



on chiral NN and 3N forces. 

The 51 < 52 Ca isotopes were produced at TRIUMF's Iso- 
tope Separator and Accelerator (ISAC) facility by bom- 
barding a high-power tantalum target with a 70 fiA pro- 
ton beam of 480 MeV energy, and a resonant laser ion- 
ization scheme [17] was used to enhance the ionization of 
the beam. 51 K was produced with 1.9 /xA of protons on a 
UC X target, and ionized using a surface source. The con- 
tinuous ion beam was extracted from the target, mass 
separated with the ISAC dipole separator system, and 
delivered to the TITAN facility. TITAN uses the time- 
of-flight ion cyclotron-resonance (TOF-ICR) [TS] method 
on singly-charged [T5H2"T] and highly-charged [22] ions 
for precision mass measurements and is capable of accu- 
rate measurements in the parts-per-billion (ppb) preci- 
sion range [23 . The system consists of a helium-buffer- 
gas-fillcd radio-frequency quadrupole |24j for cooling and 
bunching, followed by, in this case bypassed, an electron 
beam ion trap charge breeder (EBIT) [25], and a preci- 
sion Penning trap |15j . where the mass was determined. 

In the precision Penning trap a homogenous 3.7 T mag- 
netic field radially confines the injected ions, while an 
electric quadrupole field provides axial confinement. In 
order to determine an ion's mass m, the cyclotron fre- 
quency v c = qB/(2irm), where q is the charge of the ion 
and B is a homogeneous magnetic field, is determined 
from the minimum of the TOF-ICR measurement. Typ- 



ical TOF-ICR resonances for 51 K, 51 Ca, and 52 Ca are 
shown in Fig. [T] To calibrate the magnetic field, mea- 
surements with a reference ion of well-known mass were 
taken before and after the frequency measurement of the 
ion of interest. To eliminate magnetic field fluctuations 
a linear interpolation of the reference frequency to the 
center time of the measurement of v c is performed, and 
a ratio of the frequencies R — v c ,ref/vc is taken. The 
atomic mass m a of interest can then be extracted from 
the average frequency ratio R, m a — R(m a re f — m e )+m e 
where m e is the mass of the electron. 

To reduce systematic effects, such as those arising from 
field misalignments, incomplete trap compensation, etc., 
a reference with a similar mass was chosen |15j . For the 
measurements where a mass doublet was not formed a 
mass measurement of 41 K relative to 58 Ni was performed 
to investigate the mass-dependent shift. The mass of 41 K 
was in agreement with Ref. |26j . constraining the mass- 
dependent shift to be below 3.5 ppb in the frequency ra- 
tio. We include this shift as a systematic uncertainty. 
Moreover, to exclude potential effects stemming from 
ion-ion interactions from simultaneously stored isobars, 
contaminants were eliminated from the trap by applying 
a dipolar field at the mass-dependent reduced cyclotron 
frequency 1/4. for 20 ms prior to the quadrupole excita- 
tion. Dipole cleaning pushes the contaminants far from 
the precision confinement volume, greatly reducing any 
potential shifts due to their Coulomb interaction with the 
ion of interest. A quadrupole excitation time of ~ 80 ms 
was used for each of the mass measurements. In addition, 
a count class analysis |27j was performed when the count 
rate was high enough to permit such an analysis. In cases 
where the rate was too low, the frequency ratio was de- 
termined twice: Once where only one ion was detected, 
and a second time where all detected ions were included. 
The ion-ion interaction systematic uncertainty was taken 
to be the difference between these two methods and was 
100 ppb for 51 K and 750 ppb for 52 Ca. 

In previous works, mass measurements for neutron- 
rich potassium and calcium isotopes were reported up 
to 50 K and 52 Ca [25] [25], but with large uncertainties. 
While the half-lives in the vicinity of N = 32 are gen- 
erally long (tx/2 > 50 ms), the production yields in this 
high N/Z region have limited precision mass measure- 
ments until now. The mass of 51 Ca, as tabulated in the 
atomic mass evaluation ( AME2003) [25] , depends largely 
on three-neutron-transfer reactions from beams of 1 C or 
18 to a 48 Ca target 30 53]. Measurements bydifferent 
groups using the same reaction lead to differing results, 
calling into question the derived mass. Further prob- 
lems with the reaction methods include potential 40 Ca 
contamination in the target, disagreement on the num- 
ber and excitation energies of observed states, and low 
statistics. In addition, two TOF measurements of the 
51 Ca mass [351131] are in agreement with each other, but 
disagree with the values derived from the three-neutron- 
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TABLE I. Measured average frequency ratios R and mass excess (ME) with total error of 51 ' 52 Ca and 51 K, compared to the 
AME2003 mass excess. The * indicates that the AME2003 value for 51 K is based on an extrapolation. 



Isotope 


Reference 


R = v c>re f/v c 


ME (keV) 


MEameos (keV) 


ME-MEameos (keV) 


bi Ca 


b8 Ni 


0.87961718(42) 


-36338.9(22.7) 


-35863.3(93.8) 


-475.6(96.5) 


52 Ca 


58 Ni 


0.89691649(187) 


-34260.0(101.0) 


-32509.1(698.6) 


-1751.0 


52 Ca 


52 Cr 


1.00043782(158) 


-34235.8(76.4) 


same 


-1726.7 


52 Ca Average 






-34244.6(61.0) 


same 


-1735.5(701.3) 


51 K 


bly 


1.00062561(28) 


-22516.3(13.1) 


-22002.0(503.0)* 


-514.3(503.2) 




Neutron Number N 



FIG. 2. (Color online) Two-neutron separation energy San as 
a function of neutron number N for the potassium (circles), 
calcium (square), and scandium (triangles) isotopic chains. 
Points enclosed by a circle are not based completely on ex- 
perimental values in AME2003. The symbols connected by a 
dotted line are based on the TITAN mass values (mass val- 
ues for 49 ' 50 Ca and 47 - 50 K are taken from Ref. [16]), while the 
symbols connected by solid lines are those from the AME2003. 



transfer reactions. However, the uncertainties reported 
are 2-10 times larger than those obtained from the re- 
actions. All masses tabulated in the AME2003 disagree 
with the presented measurement by more than la. More 
recently, a TOF mass measurement was completed at the 
GSI storage ring [35]. This does not agree with any of 
the previous measurements and deviates by 1.3<7 from 
the value presented in this Letter. For 52 Ca, the existing 
mass value is derived from a TOF measurement [28] and 
a /3-decay measurement to 52 Sc [35]. Neither measure- 
ment agrees with our precision mass. No experimental 
data exists for the mass of 51 K. 

Our new TITAN mass measurements for 51,52 Ca and 
51 K are presented in Table [i] The mass of 51 Ca devi- 
ates by 5cr from the AME2003 and is more bound by 
0.5 MeV. We find a similar increase in binding for 51 K 
compared to the AME2003 extrapolation. For the most 
neutron- rich 52 Ca isotope measured, the mass is more 
bound by 1.74 MeV compared to the present mass table. 
This dramatic increase in binding leads to a pronounced 
change of the derived two-neutron separation energy S^n 
in the vicinity of N — 32, as shown in Fig. [2j The re- 



sulting behavior of S211 in the potassium and calcium 
isotopic chains with increasing neutron number is signif- 
icantly flatter from N = 30 to N = 32. This also dif- 
fers from the scandium isotopes, derived from previously 
measured mass excess with large uncertainties or from 
the AME2003 extrapolation. The increased binding for 
the potassium and calcium isotopes may indicate the de- 
velopment of a significant subshell gap at N — 32, in line 
with the observed high 2 + excitation energy in 52 Ca|36|. 

Three-nucleon forces have been unambiguously estab- 
lished in light nuclei, but only recently explored in 
medium-mass nuclei [SJ H3 EH] • These advances 
have been driven by chiral effective field theory, which 
provides a systematic expansion for NN, 3N and higher- 
body forces |39| . combined with renormalization group 
methods to evolve nuclear forces to lower resolution [3D] . 

We follow Ref. [5J and calculate the two-body inter- 
actions among valence neutrons in the extended pfg^/2 
shell on top of a 40 Ca core, taking into account valence- 
core interactions in 13 major shells based on a chiral 
N 3 LO NN potential evolved to low-momentum. Chiral 
3N forces are included at N 2 LO, where the two shorter- 
range 3N couplings have been fit to the 3 H binding en- 
ergy and the He charge radius. For valence neutrons 
the dominant contribution is due to the long-range two- 
pion-exchange part of 3N forces [2] |6j . In Ref. [6J , the 
normal-ordered one- and two-body parts of 3N forces 
were included to first order, and an increased binding 
in the vicinity of N = 32 was predicted. Here, we im- 
prove the calculation by including the one- and two-body 
parts of 3N forces in 5 major shells to third order, on an 
equal footing as NN interactions. This takes into ac- 
count the effect of 3N forces between nucleons occupying 
orbits above and below the valence space. For the single- 
particle energies (SPEs) in 41 Ca, we study two cases: The 
SPEs obtained by solving the Dyson equation, where the 
self-energy is calculated consistently in many-body per- 
turbation theory (MBPT) to third order; and empirical 
(emp) SPEs where the p/-orbit energies are taken from 
Ref. [1] and the 39/2 energy is set to —1.0 MeV. 

In Fig. [3] we compare the theoretical results obtained 
from exact diagonalizations in the valence space to the 
TITAN and AME2003 values for the two-neutron sep- 
aration energy and for the neutron pairing gap calcu- 
li) 

lated from three-point binding-energy differences, A„ = 
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FIG. 3. (Color online) Two-neutron separation energy 5*2n 
(top) and pairing gap a£ 3 ' from three-point binding-energy 
differences (bottom) versus neutron number N for the calcium 
isotopes. The TITAN mass values and the AME2003 values 
are shown by the symbols as in Fig. [2] The point labelled "TI- 
TAN+AME2003" is based partly on the TITAN mass values 
and complemented by the AME2003 value for 53 Ca. Theoret- 
ical predictions are shown based on chiral NN and 3N forces 
(NN+3N) in the extended pfgg /2 valence space using empir- 
ical (emp) SPEs in 41 Ca and consistently calculated MBPT 
SPEs (MBPT). 



(-l) N [B(N + 1,Z) + B(N -1,Z)- 2B(N, Z)]/2. The 
predicted S 2n is very similar for both sets of SPEs and 
is in excellent agreement with the new TITAN mass val- 
ues. For 51 ' 52 Ca, the difference between theory and ex- 
periment is only < 200 keV, but we emphasize that it 
will be important to also study the impact of the un- 
certainties in the leading 3N forces. The behavior with 

(3) 

neutron number for A„ is also well reproduced, but the 
theoretical gaps are typically 500 keV larger. Finally, we 
note also the developments using nonempirical pairing 
functionals in this region |41j . which provide a bridge to 
global energy-density functional calculations. 

In summary, the mass of 51 K has been measured with 
the TITAN facility at TRIUMF for the first time, and 
the new precision masses of 51 ' 52 Ca show a dramatic in- 
crease in binding compared to the atomic mass evalu- 
ation. The most neutron-rich 52 Ca is more bound by 
1.74 MeV, a value similar in magnitude to the deuteron 
binding energy. An increased binding around N = 32 
was predicted recently in calculations based on chiral NN 
and 3N forces 6j. The new TITAN results lead to a 
substantial change in the evolution of nuclear masses to 
neutron-rich extremes. The significantly flatter behavior 



of the two-neutron separation energy agrees remarkably 
well with improved theoretical calculations including 3N 
forces, making neutron-rich calcium isotopes an exciting 
region to probe 3N forces and to test their predictions 
towards the neutron dripline. These developments are 
of great interest also for astrophysics, as similar changes 
in heavier nuclei would have a dramatic impact on nucle- 
osynthesis j42j , and the same 3N forces provide important 
repulsive contributions in neutron-star matter [13] . 
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